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SUMMARY 


This  report  describes  progress  In  the  first  fourteen  months  of 
research  Into  erosive  burning  of  composite  solid  propellants.  Both 
theoretical  and  experimental  studies  were  undertaken.  In  the  theoretical 
study,  the  propellant  burning  process  Is  described  In  a chemically 
reacting  turbulent  boundary  layer  formed  over  the  propellant  surface. 
Governing  equations  are  formulated  and  the  order-of-magnitude  analysis 
of  these  equations  Is  carried  out,  and  boundary  conditions  are  also 
described.  Results  of  design  calculation  for  test  rig  are  presented 
in  the  section  covering  experimental  work.  The  results  show  that 
various  operating  conditions  for  velocity,  pressure  and  streamwlse 
pressure  gradient  can  be  achieved  In  the  test  rig.  Items  of  work  to 
be  performed  in  the  overall  program  are  also  listed. 
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NOMENCLATURE 


As  Frequency  factor  In  the  Arrhenius  Law 

Afc  Propellant  burning  area 

Af  Piston  area 

Q Nozzle  discharge  coefficent 

Cp  Average  heat  capacity  of  the  reacting  gases,  ^.YhC^te 

ft*  I 

Cfa  Heat  capacity  of  the  bth  species 

cs  Heat  capacity  of  solid  propellant 

Cv  Constant  volume  specific  heat  of  gases 

D Diffusion  coefficient  in  Flck's  Law 

cIap  Ammonium  Perchlorate  particle  diameter 

Activation  energies  in  the  gas  phase  reaction  and  the 
solid  surface  decomposition 

e Internal  energy  of  gases 

£ Total  body  force  in  direction 

Body  force  of  htfc  species  in  1&  direction 

H Equalizer  chamber  height 

#1  Static  enthalpy  of  gases 

hem*  Convective  heat  transfer  coefficient  for  gases  over 
propellant  surface 

Static  enthalpy  of  igniter  gases 

Static  enthalpy  of  toft  species 

f\i  Heat  of  formation  of  hlC  species 

Total  enthalpy  of  gases  at  the  exit 

Total  enthalpy  of  gases  at  the  spillage  gap 

7^  Reference  enthalpy 
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aR  Static  enthalpy  difference  across  the  boundary  layer 
Lt'\ it  Indicies  representing  coordinate  directions  1,  2 or  3 
k Species  index 

k Turbulent  kinetic  energy,  £(i4ua) 

Specific  reaction  rate  constant,  defined  in  Eq.  (23) 

L Characteristic  length  in  x-direction 

Lt  Levis  number, 

M Mass  in  control  volume 

Propellant  mass  burning  rate 
Exit  mass  flow  rate 
Igniter  mass  flow  rate 
rr>p  Piston  mass 

Mass  flow  rate  through  spillage  gap 
N Total  number  of  species 

Pressure 

Pressure  in  feed  mechanism  chamber 
R Prandtl  number, 

Ql  Latent  heat  of  vaporization  for  the  propellant 
Qs  Heat  of  reaction  at  the  propellant  surface 
W R4t*  of  heat  loss 

R Gas  constant 

% Reynolds  number  based  on  L 

Ru.  Universal  gas  constant 

't*  Burning  rate  of  solid  propellant 

Sc  Schmidt  number, 

T Gas  temperature 

If  Flame  temperature 
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T°  Standard  reference  temperature 
'Hfn  Propellant  Ignition  temperature 
Propellant  temperature 
Propellant  Initial  temperature 
Tfu  Propellant  surface  temperature 
t Time 

U/t  Reference  velocity  of  gas  flow 
H Gas  velocity  in  x-direction 

'l*  Piston  velocity 

V"  Gas  velocity  in  y-direction 

V Volume  occupied  by  gases 

\^.  Diffusion  velocity  of  tett,  species  in 
1 direction 

W Average  molecular  weight  of  gases 

Wk  Molecular  weight  of  htfc  species 

Coordinates  along  and  normal  to  propellant  surface 

Distance  moved  by  piston 

Yfe  Mass  fraction  of  htV>  species 

Change  in  mass  fraction  of  hlh  species  across  the 
boundary  layer 

Greek  Symbols 

otj  Propellant  thermal  diffusivity 
6 Boundary  layer  thickness 

Cp  Oxidizer-to-fuel  mass  ratio  in  solid  propellant 

Y Ratio  of  specific  heats 

A Thermal  conductivity  of  the  gas 

\ Thermal  conductivity  of  the  solid  propellant 
M Viscosity  of  combustion  gases 
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Stoichiometric  coefficients  for  reactants  and  products 
f Gas  density 

Reference  gas  density 

AT  Density  change  across  the  boundary  layer 
Propellant  density 

Rate  change  of  species  h due  to  chemical  reactions 
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I.  INTRODUCTION 


Thrust  and  range  are  two  important  design  parameters  of  solid 
propellant  rocket  motors.  Range  is  determined  by  burning  time  of  the 
propellant,  and  thrust  as  well  as  burning  time  depend  on  the  burning 
rate  of  the  propellant.  In  solid-propellant  rocket  motors,  it  is 
generally  understood  that  burning  rate  is  affected  by  a high-velocity  hot 
combustion  gas  flow  parallel  to  the  propellant  surface.  This  effect  is 
known  as  erosive  burning.  Thus,  ability  to  predict  erosive  burning 
effect  i8  of  prime  importance  in  the  motor  design. 

Propellant  combustion  under  non-erosive  burning  conditions  involves 
a large  number  of  chemical  and  physical  processes  including  change  of 
phase,  energy  transfer,  and  mass  flow.  Heat  transfer  from  the  hot 
combustion  gases  to  the  propellant  surface  causes  the  solid  propellant  to 
vaporize  into  gaseous  reactants.  Under  erosive  burning  conditions,  the 
high  velocity  gases  form  a turbulent  boundary  layer  over  the  propellant 
surface.  Within  the  boundary  layer,  complex  transport  processes  change 
the  temperature  and  species  distribution,  and  affect  the  chemical  reaction 
kinetics  of  gaseous  reaction  zones.  As  the  heat  transfer  rate  to  the 
propellant  surface  increases,  more  propellant  vaporizes,  and  the  gaseous 
mass  flow  rate  rises,  counter  to  direction  of  heat  transfer.  Processes 
of  mass  transfer  due  to  blowing  and  convective  heat  transfer  compete  with 
each  other  and,  depending  on  the  combustion  gas  velocity,  one  dominates 
the  other.  Consequently,  burning  rate  of  a propellant  in  erosive  burning 
conditions  can  either  increase  (more  common)  or  decrease. 

Erosive  burning  is  an  ever-present  problem  in  the  design  of  high 
thrust,  short-burning-time  solid-propellant  motors.  Current  rocket 


booster  designers  sacrifice  higher  loading  fraction  (ratio  of  propellant 
weight  to  combustion  chamber  volume),  necessary  for  better  propulsion 
performance,  because  their  knowledge  of  erosive  burning  is  Incomplete. 
Attempts  to  achieve  higher  loading  fraction  often  cause  a high  pressure 
peak  and  motor  failure  immediately  after  ignition  due  to  the  erosive 
burning  of  propellant  charge.  If  the  erosive  burning  is  serious  enough 
and  is  not  factored  into  motor  design,  a common  consequence  is  unequal 
propellant  web  burnout  with  extended  heat  exposure  of  the  chamber  wall. 

This  can  also  lead  to  chamber  failure. 

The  current  research  program  was  undertaken  to  investigate  theoretically 
and  experimentally  the  erosive  burning  problem.  Objectives  of  the 
research  program  are: 

1.  To  advance  the  state-of-the-art  of  erosive  burning  processes 
by  formulating  a theoretical  model  describing  the  important 
physical  and  chemical  phenomena  in  gaseous  and  solid  phases. 

2.  To  solve  the  theoretical  model  by  selecting  a stable,  fast 
convergent,  numerical  scheme  so  that  the  effects  of  gas  velocity, 
pressure  gradient,  and  propellant  physicochemical  characteristics 
can  be  studied. 

3.  To  verify  the  theoretical  results  according  to  the  observed 
erosive  burning  phenomena  under  various  experimental  conditions 
for  chamber  pressure,  gas  velocity,  pressure  gradients,  propellant 
ingredients,  etc. 

4.  To  develop  an  erosive  burning  rate  expression  which  will  incorporate 
all  the  important  parameters  known  to  be  affecting  the  erosive 


burning  phenomenon. 


This  report  has  been  written  to  document  progress  toward  solution 
of  the  erosive  burning  problem.  Section  II  covers  some  of  the  background 
on  erosive  burning,  and  section  III  describes  progress  in  theoretical 
work  leading  to  development  of  governing  equations,  gas  phase  kinetics, 
turbulence  modeling,  and  boundary  conditions.  In  Section  IV,  experimental 
work  related  to  design  and  fabrication  of  the  test  rig  and  instrumentation 


are  discussed 


I 

I 

II.  BACKGROUND 

Erosive  burning  behavior  has  been  investigated  by  various  methods 
as  shown  in  a comprehensive  literature  survey1  on  the  erosive  burning. 

Some  important  phenomena  observed  experimentally  and  effect  of  various 
parameters  on  erosive  burning  are  summarized  below. 

Threshold  Velocity:  The  burning  rate  of  certain  solid  propellants 

increases  only  when  the  velocity  of  combustion  gases  (or  mass 
velocity)  is  greater  than  a threshold  velocity?-7  The  threshold 
velocity  has  bean  observed  to  be  a function  of  pressure. 4 * 5 ’ 7 

Negative  Erosion:  At  low  gas  velocities  (below  the  threshold  value) 

many  investigators  have  reported  a decrease  in  the  burning  rate 
of  a solid  propellant  compared  to  strand  burning  rate. 3-8 

Effect  of  Propellant  Initial  Temperature:  An  Increase  in  the  erosive 

burning  effect  due  to  the  decrease  in  the  propellant  initial 
temperature  has  been  observed.4 

Effect  of  Gas  Temperature:  At  constant  gas  velocity  and  pressure, 

the  erosive  burning  rate  is  independent  of  the  gas  temperature. 3 

Effect  of  Pressure:  For  a given  mass  velocity,  the  erosive  burning 

is  more  pronounced  under  lower  pressures.5’6  The  influence  of 
type  of  fuel  and  mixture  ratio  on  the  erosive  burning  diminishes 
for  high  pressures.3 

Effect  of  Compressibility  of  the  Caseous  Stream:  Erosive  burning 

depends  upon  the  Mach  number  of  the  convective  stream.6 

Effect  of  Strand  Burning  Rate:  Low  burning  rate  propellants 

experience  larger  erosive  burning  effect  than  high  burning  rate 
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propellants,  and  vice  versa.2  For  high  burning  rate  propellants, 
prounounced  threshold  velocity  and  strong  influence  of  gas  mass 
velocity  on  the  erosive  burning  rate  have  been  observed.3 

Effect  of  Propellant  Characteristics:  Erosive  burning  depends  on  the 

types  of  fuel  binder  and  the  mixture  ratio  in  a propellant . 3 * 4 * 1 0 

Effect  of  Metal  Addition:  Addition  of  aluminum  has  very  little 

effect  on  the  erosive  burning.10*11 

Effect  of  Rocket  Motor  Length-to-Diameter  Ratio:  For  CMDB  Propellants, 

high  pressure  peaks  at  the  head  end  of  a rocket  motor  have  been 
observed  for  a large  length-to-diameter  ratio  grain.9 

Existing  theories  on  erosive  burning  can  be  divided  into  three  distinct 
classes  depending  upon  whether  the  theory  is  based  on  (a)  a phenomenological 
heat  transfer  theory  (without  taking  chemical  reaction  and/or  flame  structure 
into  account),  (b)  a flame  theory  (a  description  of  combustion  mechanism) 
or  (c)  an  aerothermochemlcal  analysis  including  the  considerations  of 
heat,  mass,  and  momentum  transfer  in  a chemically  reacting  boundary  layer. 

Among  phenomenological  heat  transfer  theories,  the  Lenoir-Robillard's 
theory12  is  most  popular.  The  essence  of  this  theory  is  that  erosive  burning 
rate  is  proportional  to  the  forced  convection  heat-transfer  coefficient. 

The  theory  is  analytical  to  the  extent  that  use  is  made  of  heat-transfer 
correlation  for  transpiration-cooled  surface,  and  empirical  in  the  sense 
that  the  heat  transfer  correlation  is  itself  empirical.  Compared  to 
experimental  data,  the  theory  supplies  satisfactory  agreement,  which  may 
not,  however,  arise  from  theory  validity  but  from  agreement  obtained  by 
matching  the  burning  rate  equation  to  the  experimental  data.  This  match 
is  necessary  to  determine  empirical  coefficients.  Utility  of  the  Lenolr- 
Robillard  theory  is  limited  because  rv  empirical  coefficients  must  be 
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determined  experimentally  for  every  propellant.  In  addition,  it  does  not 
include  the  effect  of  combustion  and  flame  structure. 


The  representative  work  on  erosive  burning  based  upon  flame  theory 
has  been  reported  by  Vandenkerckhove. 1 3 He  considered  the  flame  structure 
and  the  mechanism  of  solid  phase  decomposition.  However,  the  erosive 
burning  models  utilizing  combustion  mechanism  or  flame  theory  concepts 
are  handicapped  by  their  application  to  special  propellants,  sometimes, 
under  special  flow  conditions. 

In  erosive  burning  problems,  most  of  the  gaseous  reactions  occur 
within  a turbulent  boundary  layer  developed  over  the  propellant  surface 
by  the  flow  of  combustion  gases.  Therefore,  the  most  realistic  analysis 
of  the  erosive  effect  is  the  third  approach  which  employs  a description 
of  the  burning  process  (e.g.,  flame  zone  structure  and  its  interaction 
with  the  flow  field).  Early  original  work  of  Corner14  used  Prandtl- 
Karman  boundary  layer  theory  to  describe  the  flow  field.  However,  this 
work  gave  no  consideration  to  the  effect  of  mass  injection,  and  also 
assumed  the  rate-controlling  step  in  the  solid-propellant  combustion  to 
be  associated  with  gas-phase  reaction.  Modern  thought  asserts  that  such 
a rate-controlling  step  is  directly  affected  by  the  decomposition  pyrolysis 
of  a solid  propellant.  Tsuji15  analyzed  the  erosive  burning  of  a solid 
propellant  by  the  boundary-layer  approximation  in  aerothermochemistry. 

He  assumed  a laminar  boundary  layer  in  his  model  but  in  an  actual  erosive 
burning  situation  the  boundary  layer  is  turbulent.  His  model  assumed  a 
rate-controlling  step  as  the  gas-phase  chemical  reaction — contrary  to 
that  of  the  conventional  explanation  of  surface  pyrolysis.  A recent  study 
by  Razdan16  on  the  erosive  burning  of  double  base  solid  propellants  also 
considered  a chemically  reacting  boundary  layer  approach.  He  also 
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assumed  laminar  flow.  However,  Razdan's  work  takes  into  account  the 
variation  of  surface  temperature  along  the  propellant  grain  and  its  effect 
on  burning  rate.  Both  Tsujl's  and  Razdan's  models  do  not  explain  the 
negative  erosive  burning  experienced  in  some  propellants  at  low  velocities. 

In  general,  although  several  previous  theoretical  models  give  some 
insight  into  the  mechanism  of  erosive  burning,  they  Inadequately  describe 
many  of  the  phenomena  observed  experimentally.  For  example,  burning  rates 
measured  by  Marklund  and  Lake3  clearly  show  that  erosive  burning  is  affected 


by  many  parameters  such  as  oxidizer-fuel  ratio,  oxidizer  particle  size, 
pressure,  etc.  So  far,  no  work  has  been  reported  that  incorporates  the 


effect  on  erosive  burning  of  all  such  important  parameters.  Especially, 
no  work  has  been  done  to  Incorporate  a sophisticated  flame  theory  with 
turbulent  flow  processes  adjacent  to  the  burning  surface. 

In  the  current  research  program,  the  third  approach  mentioned  above 
is  being  pursued.  The  most  commonly  used  Ammonium  Perchlorate  (AP)  composite 
solid  propellants  will  be  considered  in  the  model.  In  order  to  Incorporate 
the  concepts  of  the  burning  process  and  the  flame  structure  of  the  AP  solid 
propellants  in  our  theoretical  model,  it  is  very  useful  to  review  some 
of  the  non-eroslve  (strand)  burning  mechanisms. 

Of  many  combustion  models,  three  stand  out  as  being  comprehensive 
enough  to  lead  toward  a better  understanding  of  the  strand  burning  mechanism 
of  AP  solid  propellants.  They  are  the  models  of  Summerfleld, 1 7 18  Hermance , 1 9 
and  Beckstead,  Derr  and  Price.20 


Summerf ield's  Granular  Diffusion  Flame  (GDF)  Model: 

The  model  was  originally  formulated  by  Summerfleld.17  Essential 
features  of  the  two-stage  GDF  model  are  shown  in  Fig.  1.  The  model  is  one- 
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Fig.  1 TWO-STAGE  GDF  MODEL  FOR  AMMONIUM  PERCHLORATE  TYPE  COMPOSITE  SOLID 
PROPELLANTS  (Sumraerf ield  group).18 


THREE  STEP  COMPOSITE  PROPELLANT  BURNING: 

1.  Endothermic  Zeroth  Order  Pyrolysis  of  Solid  Fuel  and  Dissociative  Subli- 
mation of  AP  to  Ammonia  and  Perchloric  Acid. 

2.  Thin  Premixed  Exothermic  NH3  + HCIO4  Reaction  Zone  I,  Assumed  to  Occur 
at  the  Surface  for  Normal  Rocket  Pressures.  Pyrolysed  Fuel  Gases  Are 
Considered  Dispersed  But  Still  Unmixed. 

3.  Exothermic  Fuel-Oxidant  Reaction,  Zone  II.  Gaseous  Fuel  Pockets  Are  Consumed 
in  Atmosphere  of  A/PA-Zone  Combustion  Products  (Dif fuslonally  and/or 
Chemical  Reaction  Controlled) 

BURNING  RATE  EXPRESSION:  1/r^  - a/p  + b/pl/3 
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dimensional  and  assumes  Chat  Che  decomposition  process  at  the  propellant 
surface  is  controlled  by  the  conductive  heat  feedback  from  two-stage 
flame  occurring  in  the  gas  phase.  Since  composite  propellant  is  of 
heterogeneous  nature,  the  fuel  and  oxidizer  gases  emerge  unmixed  from 
the  propellant  surface.  The  model  presumes  that  fuel  enters  the  flame 
zone  in  tiny  gas  pockets.  Two  reaction  zones,  ammonia/perchloric  acid 
(A/PA)  and  Fuel-oxidant,  are  assumed  to  be  parallel  to  the  propellant 
surface.  Overall  reaction  rate  of  the  fuel-oxidizer  reaction  is  determined 
by  rates  of  dlffusional  mixing  (predominant  at  high  pressures)  and  chemical 
reactions  (predominant  at  low  pressures).  It  is  quite  important  to  note 
Chat  the  pressure  dependence  of  the  heat  feedback  to  the  propellant  surface 
is  strongly  influenced  by  the  kinetics  of  the  granular  diffusion  flame 
at  pressures  above  10  atm  (even  though  the  heat  feedback  from  this  source 
is  considerably  less  than  the  heat  contribution  from  A/PA  reactions). 

The  propellant  burning  rate  relationship  for  this  model  is  based  on  inte- 
grated, steady-state,  one-dimensional  energy  equations  for  the  various 
stations  in  the  flame  zone.  The  proposed  burning  rate  equation  is  a 
function  of  pressure,  chemical- react ion  time  parameter  and  diffusion  time 
parameter  (which  in  turn  is  a function  of  oxidizer  particle  size).  These 
parameters  are  experimentally  evaluated,  and  correlation  of  the  burning  rate 
equation  was  extremely  good  for  a variety  of  propellants  over  the  pressure 
range  of  1 to  100  atmospheres. 

Hermance’s  Heterogeneous  Reaction  Model: 

The  Hermance  model  is  based  on  a detailed  combination  of  the  steady- 
state  decomposition  processes  of  the  fuel  and  oxidizer.  Three  regions  of 
propellant  combustion  are  identified:  propellant  surface,  the  gas  phase 

flame  zone,  and  the  region  between  the  surface  and  the  flame.  The  surface 
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chemical  processes  are  endothermic  fuel  pyrolysis,  exothermic  oxidizer 
decomposition  and  an  exothermic  heterogeneous  chemical  reaction  between 
fuel  binder  and  decomposed  oxidizer  in  small  regions  surrounding  individual 
oxidizer  crystals.  Each  of  these  reactions  produces  a mass  flux  from  the 
propellant  surface.  The  area  on  which  the  surface  reaction  occurs  is 
calculated  by  postulating  that  an  oxidizer  crystal  decreases  in  size 
during  decomposition.  This  size  reduction  produces  a shallow  fissure 
between  the  oxidizer  crystal  and  the  fuel  binder.  This  fissure  is  the 
region  where  the  heterogeneous  reaction  takes  place.  By  using  geometric 
and  statistical  methods,  the  size  reduction  of  the  oxidizer  particles  is 
expressed  as  a function  of  oxidizer  volume  fraction,  burning  rate,  pressure, 
and  oxidizer  particle  distribution.  Although  model  is  more  detailed  than 
GDF  model,  the  predicted  burning  rates  from  this  model  were  higher 
compared  to  experimental  data.  Further,  the  experimental  evidence  of  Boggs34 
argues  against  the  fissure  postulation  in  this  model. 

The  Beckstead,  Derr,  Price  (BDP)  Multiple  Flame  Model: 

In  this  model  the  existence  of  three  flames  in  the  region  of  each 
exposed  oxidizer  particle  have  been  postulated.  Essential  features  of 
the  BDP  model  are  shown  in  Fig.  2.  Unlike  GDF  model,  AP  monopropellant 
flame  is  considered  distended  from  the  propellant  surface.  Three  flames 
considered  are:  (1)  a primary  flame  between  the  decomposed  products  of 

the  binder  and  the  oxidizer,  (2)  a premixed  oxidizer  monopropellant  flame, 
and  (3)  a final  diffusion  flame  between  the  products  of  the  oxidizer 
monopropellant  flame  and  the  decomposed  products  of  the  binder.  Energy 
feedback  to  the  propellant  surface  arriving  from  various  flames  is  considered 
separately  and  then  superimposed.  For  high  pressures,  the  reaction  path 
A in  Fig.  2 is  considered  more  important  than  the  reaction  Path  B.  The 
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Fig.  2 POSUTLATED  FLAME  STRUCTURE  FOR  AN  AP  COMPOSITE 
PROPELLANT  (BDP  Model)20 


TWO  PROPOSED  REACTION  PATHS: 

Path  A: 

I.  AP  Monopropellant  Flame:  NH3  + HCIO4  Reaction.  Independent  of  mixing 

but  kinetics  dependent.  This  flame  is  not  considered  as  occuring  at  the 
surface  for  low  pressures. 

II.  Final  Diffusion  Flame:  Oxygen  rich  products  of  reaction  I react  with 

the  fuel-rich  binder  pyrolysis  products.  Dependent  on  mixing  only. 

Path  B: 

III.  Primary  Flame:  Fuel  Pyrolysis  proudcts  react  with  HCIO4  products  giving 

combustion  products.  Dependent  on  mixing  and  kinetics.  Unimportant 

at  high  pressures. 


a 
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physical  nature  of  the  oxidizer  particle  size  and  the  burning  surface  has 
been  described  statistically,  and  the  modal  Is  capable  of  describing 
various  qualitative  characteristics  observed  In  composite  propellants 
(e.g.,  plateau  burning). 

None  of  the  solid  propellant  flame  models  discussed  above  will  predict 
the  burning  rate  of  a composite  propellant  without  prior  knowledge  of 
propellant  combustion  characteristics.  In  each  case,  various  parameters 
in  the  model  which  are  not  accurately  known  must  be  varied  to  match  the 
burning  rate  characteristics  with  the  experimental  results.  Nevertheless, 
for  the  present  problem,  these  models  suggest  a strong  interaction  of 
the  flame  zone  structure  with  the  flow  field.  For  example,  a preliminary 
calculation  based  on  experimental  data  of  references  3,  18,  21  was  made 
and  results  plotted  in  Fig.  3.  In  these  calculations,  it  has  been  assumed 
that  a viscous  sublayer  exists  and  that  the  gas  flow  does  not,  significantly, 
change  the  diffusion  flame  thickness.  This  plot  shows  that  the  diffusion 
flame  is  outside  the  viscous  sublayer  for  a high  velocity  combustion  gas 
stream.  It  is  quite  evident  that  the  role  of  turbulence  in  diffusional 
mixing  and  in  rate  of  heat  transfer  to  the  propellant  surface  will  be 
quite  important  in  erosive  burning. 

The  surface  roughness  is  another  factor  which  should  be  taken  into 
account  in  erosive  burning.  At  normal  rocket  pressures  the  roughness  could 
be  of  the  order  of  10  micron.22  This  again  depends  on  propellant  formulation 
and  in  many  propellants  the  surface  roughness  can  be  as  large  as  the 
oxidizer  particle  size.  Therefore,  the  Reynolds  number  based  on  the 
roughness  size  may  play  an  important  role  in  turbulence  characteristics 
near  the  propellant  surface. 
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Height,  pm 


60 


Pressure,  atm  — *- 

Fig.  3 Comparison  of  Viscous  Sublayer  Thickness  and  Diffusion 
Flame  Thickness  (viscous  sublayer  calculation  based 
on  experimental  data  on  erosive  burning  of  Ref.  3, 
and  the  data  on  friction  coefficient  of  Ref.  21) 
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III.  THEORETICAL  MODELING 

The  burning  rate  of  a propellant  primarily  depends  on  the  rate 
of  heat  feedback  from  reacting  gases  to  propellant  surface.  Heat  feedback 
is,  in  turn,  controlled  by  various  mechanisms  which  depend  on  parameters 
such  as  turbulence  level  in  the  flow,  flame  height,  viscous-sublayer 
thickness,  pressure,  etc.  A flow  chart  (Fig.  4)  is  constructed  to  illustrate 
our  understanding  of  the  erosive  burning  mechanism.  The  flow  chart  is 
based  on  previous  studies  of  many  investigators . 1_ 2 0 ’ 3 6 Although  inconclusive, 
the  flow  chart  is  useful  in  the  development  of  a comprehensive  theoretical 
model.  The  influence  of  the  above  mentioned  parameters  can  be  successfully 
studied  by  considering  heat  feedback  to  the  propellant  surface  within  a 
turbulent  boundary  layer. 

Flow  of  high  velocity  combustion  gases  over  a rocket  motor  solid- 
propellant  surface  forms  a turbulent  boundary  layer.  The  essential 
features  of  a rocket  motor  are  shown  in  Fig.  5.  Most  reactions  occur 
within  the  boundary  layer,  and  the  flow  field  interacts  strongly  with 
the  flame  zone  structure.  The  third  approach  mentioned  earlier  (i.e., 
an  aerothermochemical  analysis  in  chemically  reacting  turbulent  boundary 
layer,  employing  a description  of  the  burning  process)  is  therefore  the 
most  logical  analysis  for  the  erosive  burning  problem.  In  the  theoretical 
model,  a steady,  two-dimensional  compressible,  turbulent  boundary  layer 
over  a flat  plate  has  been  considered  (Fig.  6). 

III. A.  Governing  Equations 

To  formulate  the  theoretical  model,  we  start  with  the  general 
conservation  equations*  for  a reacting  compressible  fluid  flow. 

*See  nomenclature  for  the  definition  of  various  mathematical  symbols. 
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TUB  FLOW  OB  HIGH  VELOCITY, 
PRODUCT  CAS  OVEt  MJRMING 
PROPELLANT  SURFACE  PRODUCES 
A TURBULENT  BOUNDARY  LAYER. 


FIG.  4 PLOW  CHART  FOR  EROSIVE  BURNING  MECHANISM  OF  COMPOSITE  SOLID  PROPELLANTS 
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Fig.  5 Schematic  Diagram  of  a Solid-Propellant  Rocket  Motor 
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Fig.  6 Reactive  Turbulent  Boundary  Layer 
Formation  on  the  Test  Plate 


Mass  conservation: 


K +Cfli<),i=o  (1) 

Momentum  conservation: 

^(fu.)  + (Pa.uj)t.  = -*fi  + Tijfi  + f*  (2) 

Species  conservation: 

+ (fmfc  + (PYhVki)>i=  , f»=u. n (3) 

Energy  conservation: 


- $ - M*  = - PWi^PYOA} 

ft®! 


+ Tij  U’irj 

(4) 

where : iv 

,T 

f<>ZY.k  - 

K = + 1 ^ 

(5) 

I 

T* 

£ 

II 

> 

p 

(U. . 

u“5  - 

(6) 

Equation  of  state: 

p = fRT 

(7) 

In  these  equations , each  variable  is  replaced  by  its  mean  and  fluctuating 

part  (e.g.,  Ui=Ui  + U/i  , fi  - J,  + ^ , etc.)*  While  doing  so,  the  diffusion 

velocity,  \^_  , appearing  in  equation  (3)  and  (4)  is  replaced  by  using 

Fick's  Law  of  Binary  Diffusion,  viz.  ^ . After  neglecting 

'te 

the  body  force  terms,  each  equation  is  averaged  and  the  following  governing 
equations  for  the  mean  quantities  are  derived: 

Continuity  Equation: 

$[+  + = 0 (8) 
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Momentum  Equation: 

¥ gvi  + .(?*;  4-  (fttj  4 fuj)  = - it  + (TLtj  - cfu/5ui).  O) 

Species  Equation: 

f f'  * !/?  ♦ = (f«V  - * Ow  do 

Energy  Equation: 

ff  + # * + 7£)l,i-gt  - a, 

= — (pttiili  + £4  + ( u - J . 

-4  Tii  l Ji,j  - FK.i  + xtj  Ui,f  (11) 

Reynolds  Stress  Equation: 

This  equation  is  derived  as  follows:  “Subtract  momentum  equation 
for  mean  quantities  from  the  momentum  equation  for  Instantaneous 
quantities  and  multiply  by  (1< . To  the  resulting  equation,  add 
a similar  equation  with  lndicies  1 and  l interchanged.  Then 
average  the  whole  equation,  and: 

+ ypKuj)  + (ruj4f^)(uTuJ^ 

-4  (fUj/ui  u(fj  -4  (fu^'ul  Ui,j  + (fa'iUia'j),- 

+ ( Oj  C^uTT*  uMp'iij  - 

= ~ ^ ^0,)  "4  u*)»-i3  ■*"  *P  Cdi»£  + H?,/)  4 ( 4 (12) 


Boundary  Layer  Equations  and  Order  of  Magnitude  Analysis: 

Following  assumptions  are  Introduced  into  Eqs.  (8)  to  (12): 
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1.  Average  quantities  are  steady  state 

2.  Mean  flow  Is  two  dimensional 

3.  For  the  Mach  number  range  ^ 1 


^a_«  . 


- bH 


and 


f'a'V-' 

p UFF 


« 1 


4. 

5. 


According  to  Laufer,2**  who  based  his  conclusion  on  the  experi- 
mental results  of  Kiatler,22  the  density  fluctuations  have  a kinematic 
rather  than  a dynamic  effect  on  the  turbulence. 

There  is  no  reaction  generated  turbulence25 

T ~~  0(i) 


6.  Le  = 1 

7.  | ( Tij  UiUif'  + UiUiP'uJ-  - C5u5  f^i)  <^<  Jjjs? 

°i  * 

8.  Pe  » i » ’ 4«1 


Following  the  regular  procedure  for  order  of  magnitude  analysis  (see 
Appendix  A),  dominant  terms  in  the  Eqs.  (7)  - (12)  are  retained.  Turbulent 
kinetic  energy  equation  (4=0  and  Reynolds  shear  stress  equation  ( <-L  t ^=2) 
are  obtained  from  Eq.  (12).  Finally,  the  two  dimensional  steady-state 
turbulent  boundary  layer  equations  for  a reactive  turbulent  compressible 
flow  are  obtained: 


Continuity  Equation: 

fc(W>  + O 

where 

3 ??+fV 


(13) 


Momentum  Equation: 


Species  Equation: 

pag  4-  - (W£j  + c\ 

Energy  Equation: 

fn|6  + fjF»|  = 


(15) 


(16) 


Contribution  of  , ^(4^  and  U'i,j  in  Eq.  (16)  is 

generally  believed  to  be  small  compared  with  the  convection  of  £ and  is 
therefore  neglected  in  the  subsequent  analysis. 

Turbulent  kinetic  energy  equation: 


' a# 


- r^F(k'  ■"  f)vrj  -Z4 


(17) 


Reynolds’  shear  stress  equation: 


+ *>'(?-}'  * 


—2ju 


i>u'  ay-' 

*1  ^3 


(18) 


Equation  of  State: 


f = ^(F*+P'A') 


(19) 
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Heat  Conduction  Equation  in  the  Solid  Phase: 

It  is  assumed  that  the  heat  conduction  into  the  solid  propellant 
is  dominant  in  a direction  normal  to  the  burning  surface.  In 
a coordinate  system  attached  to  the  burning  surface,  the 
temperature  distribution  in  the  solid  propellant  at  a given  x 
location  along  the  surface  is  governed  by: 

’ where <^  °'  (20) 

III.B.  Modeling  of  Gas-Phase  Chemical  Reactions 

When  a composite  solid  propellant  burns,  solid  fuel  and  oxidizer 
particles  transform  into  gases.  The  gases  may  react  in  several  steps. 

However,  in  the  present  work,  the  following  single-step  chemical  reaction 
is  assumed: 

0 * F — p (21) 

where  0 and f represent  the  oxidizer  and  the  fuel  gases  respectively,  and  P 
represents  the  product  gases.  In  view  of  the  already  complicated  turbulent 
transport  phenomena,  the  above  assumption  is  justified. 

For  this  single-step  reaction,  production  of  Infc  species  due  to  the 
chemical  reaction  is  given  by  the  following  equation: 

= W„  C V<>  *.ff  (22) 

h»i 

where  the  specific  reaction  rate  constant,  , is  given  by  the  Arrhenius 
law: 

A , (23) 

in  which  the  frequency  factor:  A = BT11  , 

Usually  in  a homogeneous  reaction  B is  a constant.  However,  for 
diffusion  controlled  gas-phase  reactions  of  a burning  solid  propellant, B 
will  not  be  a constant.  It  is  assumed  that,  at  least,  B has  the  following 
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functional  dependence: 


B=B(<*ap  »<#>)  (24) 

From  Eq.  (22),  the  rate  of  consumption  of  oxidizer  and  fuel  species  can  be 
represented  by: 

u.  - - 4 ’ C*  - - & (’Yy. 

Applying  Reynolds'  decomposition  to  (J0  (=  6J0  + oJi),and  then  averaging, 
we  have: 


= - ^ ^2?Y.  (25) 

A similar  equation  is  valid  for  6JF  with  VJF  replaced  by  Wo,  in  Eq.  (25). 

In  the  above  equation,  fluctuations  in  &s  are  assumed  negligible 
because  temperature  fluctuations  will  not  introduce  any  significant 
fluctuations  in  for  large  values  of  E-o-j.  in  Eq.  (23).  In  Eq.  (25), 
third  and  fourth  order  correlations  have  been  neglected.  It  has  been 
common  to  neglect  all  the  turbulence  correlations  on  the  light  hand  side  of 
Eq.  (25).  However,  this  is  incorrect  when  the  reaction  rates  are  fast,28 
i.e.,  when  the  chemical  reactions  are  diffusion  controlled  which  is 
valid  for  burning  composite  solid  propellants. 

Relative  importance  of  various  correlations  in  Eq.  (25)  is  currently 
being  investigated.  It  is  believed  that  the  correlation  Yo' V/  should  be 
related  to  the  fuel  pocket  size  (GDF  theory) , which  is  dominant  in  diffusion 
controlled  chemical  reactions. 

III.C.  Turbulence  Closure 

To  solve  Eqs.  (13)  - (19),  a number  of  turbulence  correlations 
(••«•»  (pvY . (PVO'fi  t «tc.)  have  to  be  modeled.  Some  authors35 
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have  developed  conservation  equations  for  such  correlations  and  thereby, 
introduced  higher  order  correlations.  The  higher  order  correlations  are 
usually  modeled  in  terms  of  the  known  variables  and  some  unknown  constants. 
These  constants,  in  general,  are  evaluated  by  matching  the  theoretical 
and  experimental  results.  However,  in  our  current  approach,  excess 
complexity  is  avoided  in  favour  of  simpler  and  faster-to-execute  formulation. 
Also  it  is  understood  that  chemical  reactions  in  the  diffusion  flame 
play  an  important  role  in  erosive  burning.  Therefore  it  would  be  more 
meaningful  to  describe  the  effects  of  turbulence  on  the  chemical  reactions 
in  the  closure  consideration,  rather  than  solve  numerous  conservation 
equations  for  higher  order  turbulence  correlations.  For  the  gas  phase, 
only  Eqs.  (13)  - (19)  are  treated  as  the  governing  equations  in  the  theoretical 
formulation.  Terms  to  be  modeled  for  turbulence  closure  are  given  below. 


1.  Velocity  diffusion  terms, 


) 

2.  Pressure  diffusion  terms,  + |-(t>Li-j) 

3.  Redistribution  terms, 

r v 3k*  / 

4.  Dissipation  terms,  /X 


5.  Correlation  appearing  with  Production  terms, 

6.  Turbulence  produced  heat  transfer,  (P<V)£ 


7Z7 


7.  Turbulence  produced  species  transfer,  (fU/)Yh 


8.  Species  correlation  term,  Y0'yF' 

Currently,  the  modeling  effort  of  the  various  terms  listed  above  is 


in  progress.  For  some  of  the  terms,  existing  models  of  Rotta29  and 


- 23  - 


jL 


Donaldson  28  will  be  examined  for  a possible  extension  to  the  present 
compressible  turbulent  boundary  layer. 


III.D.  Boundary  Conditions 

Mass  balance  at  the  propellant  surface: 

Consider  a thin  Control  Volume  (C.V.)  around  the  solid-gas  interface 
(Fig.  7),  mass  balance  for  species  k can  be  written  as 


+ 


(26) 


TMass  Rate  of 
hth  species  of 
solid  propellant 
Into  C.V. 


Mass  Rate  of  h.th 
Species  of  Gases 
Transported  out 
of  C.V.  by  surface 
blowing. 


Mass  Rate  of  feth 
Species  of  Gases 
Transported  out 
of  C.V.  by  diffusion. 


Using  (VV~)0+  1x1  E9*  (26)»  we  have 


(27) 


Yf(«»o*)  = 


\7n.o*) 


Vu.o*)  (»*o4) 


. YPc«.o*) 


(28) 


Since  no  product  species  can  diffuse  into  the  solid,  we  have 

Y,C*»cn  + = i 

or  . 

t ♦9 

Also  from  Fick's  law,  the  diffusion  velocity  is  given  by 


and 


^Z(*,o*)  - - —^2 — £!f| 


(29) 

(30) 


(31) 


io 


(32) 


Energy  balance  at  the  propellant  surface: 

Consider  a similar  solid-gas  Interface  control  volume  as  before; 
the  energy  balance  at  the  propellant  surface  (Fig.  8)  Is  given  below. 


A‘^0- 

= 

+ 

fs^bCQs-QL)|o 

Rate  of  Heat 

Rate  of  Heat 

Net  Rate  of  Heat 

Flux  Conducted 

Feedback  from 

Produced  by 

into  the 

gases  to 

Exothermic  Reaction 

Propellant . 

surface. 

inside  C.V. 

+ 

^b%.CCoYo 

<r 

Rate  of  Energy  carried  into 
the  C.V.  by  Propellant  Feeding. 


— f (v-  +VF)  Yf 

Rate  of  Energy  Transported 
away  by  Convection  and 
Diffusion  of  Fuel  Species. 

- (33) 

I o 

Rate  of  Energy  Transported 
away  by  Convection  and 
Diffusion  of  Oxidizer  Species. 

At  the  propellant  surface,  energy  carried  by  convection  of  product 
species  and  energy  carried  by  diffusion  of  product  species  balance 
each  other.  Therefore,  there  is  no  net  contribution  of  product  species 
in  the  energy  balance  equation. 

Substituting  Eqs.  (27)  and  (28)  into  Eq.  (33),  and  neglecting  the 
contribution  due  to  diffusion  ( \^<<V  ) , we  have: 

\^\  _ = * t 00  + T*u>  (cj  - CO]  (34) 
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In  addition  to  Eq.  (34),  the  other  boundary  condition  for  the  heat 


conduction  equation  in  solid  phase  is: 


Tj,  ( x , - •*  ■)  - 

(35) 

The  other  boundary 

conditions  needed  at  the  propellant  surface  are 

as 

follows : 

0Ux,o)  = o 

(36) 

Ps 

VCx.o)=  ■ b 

(37) 

P(x,o) 

where,  the  solid  propellant  burning  rate  is 

given  by  the  Arrhenius 

law: 

(38) 

(39) 

The  boundary 

conditions  on  the  adiabatic  surface  (Fig.  6)  of 

the 

leading  edge  are: 

U Cx,o*)  = o 

(40) 

V-(x,o+)  = o 

(41) 

O 

1! 

■f 

o 

^ol/o 

(42) 

>£(Xi  <0  = o 

(43) 

Yp(X,  0*)=  0 

(44) 

The  boundary  conditions  at  the  free  stream  are: 


CL  (x,  ®°")  = U-00  (45) 

T(x,-)  = Tot*)  (46) 

$(x»®0  = lj.00  (47) 

YoCx.oo-)  - o (48) 

YF  (*.•*)  = ° 
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(49) 


To  solve  the  differential  equations  for  k and  uV'  , we  also 
need  boundary  conditions  on  them.  These  boundary  conditions  are  currently 
under  investigation. 
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IV . EXPERIMENTAL  WORK 


The  objective  of  this  portion  of  the  work  is  measurement  of 
the  burning  rate  of  a solid  propellant  ant  study  of  related  variables. 
Experimental  data  obtained  can  then  be  used  to  verify  the  theoretical 
model . 

Two  requirements  were  imposed  on  design  of  the  experimental  apparatus 
Experiments  should  be  compatible  with  the  theoretical  model,  and  conditions 
of  test  propellant  burning  should  be  similar  to  that  surrounding  combustion 
of  a typical  solid  propellant  motor.  Therefore,  a propellant  sample  should 
burn  at  high  pressures  (50-100  atm) , a turbulent  boundary  layer  should  be 
formed  over  the  propellant  sample  by  high  velocity  (100-200  m/sec)  gas 
flow,  and  the  propellant  sample  surface  should  be  flat  during  combustion. 

IV. A.  Burning  Rate  Measurement 

A promising  method  for  measuring  burning  rate  is  the  Laser-Photodiode 
servome.chanism  technique.  It  has  not  been  broadly  used  in  erosive  burning 
experiments.  Advantages  of  this  technique  are:  the  burning  rates  measured 

are  instantaneous  and  accurate,  and  the  results  of  th'i  burning  rates  are 
directly  available  through  an  electronic  recording  system.  The  basic 
components  used  in  the  technique  are:  a He  - Ne  Laser,  a light-sensitive 

photodiode,  a control  amplifier,  a recording  system  and  a stepper  motor. 
Schematic  diagram  of  the  Laser-Photodiode  servomechanism  technique  is 
shown  in  Fig.  9.  As  the  surface  of  the  burning  solid  propellant  strip 
recedes,  the  laser  beam  (previously  blocked  by  the  propellant)  is  detected 
by  the  photodiode  to  generate  a small  current  signal.  This  signal  is 
amplified  in  the  control  amplifier,  and  supplied  to  the  stepper  motor, 
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which  advances  the  propellant  strip  counter  to  the  receding  propellant 
surface.  If  the  stepper  motor  advances  the  propellant  too  quickly,  the 
laser  beam  will  be  blocked,  the  signal  to  the  motor  discontinues,  and 
propellant  strip  feeding  stops.  Thus,  feeding  rate  is  controlled  by 
burning  rate.  Both  surface  recession  from  burning  and  surface  advancement 
by  the  stepper  motor  are  very  fast.  With  respect  to  a stationary  observer, 
the  surface  of  the  propellant  is  maintained  at  a fixed  location.  The 
feeding  rate  by  the  stepper  motor  is  recorded  by  the  recording  system  at 
all  times  during  the  experiment,  thus  giving  the  instantaneous  burning  rate. 

As  mencioned  earlier,  a flat  solid  propellant  surface  is  necessary 
for  validation  of  the  theory,  which  is  modeled  for  a flat  plate  turbulent 
boundary  layer. 

When  a long  solid  propellant  surface  burns  in  a gas  stream,  the 
surface  does  not  regress  equally  at  all  points  along  the  surface  in  the 
flow  direction.  Consequently,  the  surface  no  longer  remains  flat.  This 
problem  is  solved  by  designing  an  experimental  setup  which  will  feed 
four  solid  propellant  strips  at  different  rates  such  that  the  flat  surface 
is  maintained.  This  is  achieved  by  the  experimental  setup  shown 
schematically  in  Fig.  10.  Four  laser  beams  together  with  four  stepper 
motors  control  the  position  of  the  four  different  propellant  surfaces  and 
maintain  each  surface  at  the  same  level. 

IV . B . Test  Rig  Design 

The  design  of  various  components  of  the  test  rig  depends  on  the 
gas  dynamic  operating  conditions  (pressure,  temperature  and  velocity). 

Basic  equations  used  for  the  design  of  the  test  rig  are  given  in  Appendix 
B.  A parametric  study  for  the  effects  of  various  geometric  inputs  of  the 
, test  rig  on  the  solutions  of  these  equations  is  also  presented  in  Appendix  B. 

I 
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Various  materials  and  dimensions  of  the  test  rig  were  selected  on  the 
basis  of  the  parametric  study  and  the  requirements  of  high  pressure  and 
high  temperature  gas  flow.  A schematic  diagram  of  the  test  rig  is 
shown  in  Fig.  11.  The  side  (multi  cross-sectional)  view  is  shown  in  Fig.  12. 
Description  and  function  of  each  important  component  is  given  below: 

1.  Driving  Motor:  For  meaningful  experimental  results,  the  test 

conditions  for  the  burning  of  propellant  strips  should  be  made 
similar  to  the  burning  of  a propellant  in  a rocket  motor.  Conse- 
quently, a high  pressure,  high  velocity,  hot  combustion  gas  flow 
is  needed.  Such  a gas  flow  is  generated  in  the  driving  motor  by 
burning  a 12-inch-long,  star  shaped,  solid  propellant  grain.  The 
propellant  grain  (N-4)  is  the  same  as  used  in  the  5-inch  sidewinder 
1A  motor,  and  has  neutral  burning  characteristics  (constant  mass 
burning  rate)  which  help  to  maintain  a steady  state  pressure  in 
the  driving  motor  after  a short  ignition  transient  interval. 

The  solid  propellant  grain  is  held  in  a 304  stainless  steel 
cylindrical  container,  designed  to  withstand  pressures  up  to 
800  atm. 

To  Ignite  the  propellant  charge  in  the  driving  motor,  a 
pyrotechnic  igniter  (MK  125  MOD  5)  which  contains  a mixture  of 
black  powder  and  magnesium,  will  be  used.  The  igniter  is  set  off 
by  an  electric  squib. 

After  the  ignition  of  the  propellant  grain,  the  product  gases 
flow  out  of  the  driving  motor  into  the  test  chamber  through  a 
nozzle.  This  nozzle  (304  stainless  steel)  is  designed  to  converge 
from  a square  cross  section  (80  cm2)  at  the  end  of  driving  motor 
to  a rectangular  cross  section  (20  cm2)  at  the  entrance  of  the 
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test  chamber.  To  minimize  the  heat  losses  through  the  nozzle, 
the  internal  surface  of  the  nozzle  is  coated  with  a thin  insula- 
tion layer  of  zirconium  oxide  (Union  Carbide  LZ-2). 

Test  Chamber:  Several  important  elements  are  designed  to  form 

the  39  cm  long  test  chamber  having  a rectangular  cross-section 
of  7 cm  x 2 cm. 

The  elements  are  described  below. 

(a)  An  interchangeable  wedge  shaped  stainless  steel  leading 
edge  which  is  thermally  insulated  with  LZ-2  coating  to 
prevent  heat  losses.  The  leading  edge  is  attached  to  a 
stainless  steel  flat  plate  through  which  four  test  propellant 
strips  (7.5  cm  x 3 cm  x 2.5  cm)  are  fed  by  the  stepper  motors. 
The  propellant  strips  are  thermally  insulated  from  the 

flat  plate  by  using  a quartz  lining  around  the  strips. 

This  is  necessary  to  avoid  any  unequal  burning  due  to  high 
rate  of  heat  conduction  to  the  propellant  edges  in  contact 
with  plate.  The  length  of  the  leading  edge  (11  cm)  is 
selected  to  ensure  that  the  boundary  layer  is  turbulent 
over  the  propellant  surface. 

(b)  A spillage  channel  through  which  about  10  percent  of  the 
product  gases  flow  out  of  the  test  rig.  This  ensures  that 
the  boundary  layer  starts  to  develop  from  the  beginning 

of  the  leading  edge. 

(c)  An  interchangeable  top  plate  is  used  to  vary  the  channel 
height  in  the  test  section  to  change  gas  velocity.  Pressure 
gradient  can  also  be  controlled  by  using  a tapered  top 


plate. 


(d)  A convergent-divergent  interchangeable  exit  nozzle  is 
designed  to  control  the  mass  flow  and  the  gas  velocity  in 
the  test  chamber.  The  exit  nozzle  assembly  contains  a burst 
diaphragm  designed  to  rupture  at  a given  critical  pressure. 

A small  nozzle  with  a throat  area  of  about  one  tenth  of 

the  exit  nozzle  is  used  at  the  end  of  spillage  channel.  The 
small  nozzle  assembly  also  contains  a burst  diaphragm 
designed  to  rupture  at  a pressure  slightly  less  than  the 
critical  pressure  of  the  main  exit  diaphragm.  This  is 
to  ensure  the  opening  of  the  spillage  channel.  Through  the 
use  of  these  diaphragms,  pre-pressurization  of  the  test 
chamber  can  be  achieved. 

(e)  A plexiglass  window  assembly  is  made  of  a sacrificial 
plexiglass  plate  and  a main  plexiglass  window.  This  assembly 
is  installed  in  the  front  cover  of  the  test  chamber.  To 
enable  the  laser  beams  to  pass  through  the  test  section 

and  to  be  sensed  by  photodiodes,  four  small  cylindrical 
quartz  rods  are  located  in  the  back  plate  of  the  test  chamber. 

3.  Feed  Mechanism  Chamber:  The  chamber  houses  four  stepper  motors 

and  four  feed  screws.  Rotational  motion  of  the  feed  screws 
(attached  to  the  stepper  motors)  is  converted  to  linear  vertical 
motion  of  a rectangular  element  through  a threaded  nut  (see 
Fig.  13).  The  propellant  strip  is  carried  by  each  of  the  rectangular 
elements.  Light  weight  aluminum  was  chosen  as  the  material  for 
the  rectangular  element.  Four  stainless  steel  feed  screws  (3/8 
in.  fine  SAE  24  threads/in.)  are  used  to  feed  the  propellant  strips. 

The  torque  requirement  of  the  stepper  motor  depends  on:  the 
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total  load  on  the  motor,  friction  and  efficiency  of  the  feed 
screw.  To  minimize  the  torque  requirement  of  the  stepper  motors, 
the  pressure  difference  between  the  test  chamber  and  the  feed 
mechanism  chamber  should  be  kept  minimal.  This  is  achieved  by 
pressure  equalizer  chamber,  discussed  below. 

4.  Pressure  Equalizer  Chamber:  The  pressure  equalizer  chamber  consists 

of  a piston  and  a cylinder.  The  piston  is  connected  with  a spring 
and  a damper  to  reduce  oscillations  of  the  piston.  The  chamber 

is  connected  to  the  test  chamber  and  to  the  feed  mechanism 
chamber.  Higher  pressure  in  the  test  chamber  will  push  down  the 
piston  compressing  the  gas  below  it  and  increasing  the  pressure 
in  the  feed  mechanism  chamber.  Consequently,  the  pressure 
difference  between  the  test  chamber  and  the  feed  mechanism  chamber 
is  reduced. 

5.  Instrumentation: 

(a)  Stepper  motors:  A stepper  motor,  producing  150  pulses/sec 

with  an  index  angle  of  45°  gives  a linear  feeding  rate  up 

to  2 cm/sec  when  24  threads/in.  screw  is  used.  This  feeding 
rate  is  within  the  typical  values  for  erosive  burning  rates. 
Stepper  motor  control  unit  is  used  to  supply  the  necessary 
voltage  and  current  to  the  stepper  motor.  Both  these  Items 
are  obtained  from  Clifton  Precision  Company.  The  burning  rate 
of  a solid  propellant  strip  can  be  measured  directly  by 
recording  the  feeding  rate  of  the  stepper  motor. 

O 

(b)  He-Ne  laser:  A 6328  A and  1-2  mw  laser  (Spectra-Physics  model 

136)  is  used.  It  has  adjustable  power  output,  precision  mounting 
surfaces,  and  proper  dimensions  for  the  current  design. 
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(c)  Beam  splitter:  The  total  number  of  lasers  needed  is 

reduced  by  half  by  using  beam  splitters.  These  are  used 
to  split  two  original  laser  beams  into  four  beams  of  equal 
power.  Ealing  Corporation's  two  splitter  cubes  having  equal 
transmission  and  reflection  with  two  reflection  prisms  are 
used  for  this  purpose. 

(d)  Band-pass  filter:  Before  the  laser  beam  reaches  a photodiode, 

a band-pass  filter  with  a narrow  frequency  range  is  used  to 
filter  out  any  signal  produced  by  light  emission  of  combustion 
gases.  The  band-pass  filter  (Oriel  Corporation  model  G-577- 

O 

6328)  having  a band  width  of  10  A is  used  for  isolating  the 

O 

6328  A laser  beam  signal. 

(e)  Pressure  transducers:  These  are  used  to  measure  the 

pressure  at  the  upstream  and  downstream  locations  in  the 
test  chamber.  The  pressure  transducers  (Kistler  model  601  B) 
having  pressure  range  up  to  1,000  atm  have  bee  selected. 

These  transducers  are  equipped  with  water  cooled  adapters 
(Kistler  Model  628  C). 

(f)  Thermocouple  probes:  Omega  Engineering  Tungsten  - 5 percent 

Rhenium  vs.  Tungsten  - 26  percent  Rhenium  thermocouple  probes 
having  Beryllia  insulation  have  been  selected.  Temperatures 
up  to  2600°K  can  be  monitored. 

Other  instrumentation  includes  control  amplifier, 
photodiodes,  recording  system,  lead  feedthroughs,  laser 
exciters,  stepper  motor  controls,  etc. 
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During  the  reporting  period,  the  following  items  have  been  completed 
as  part  of  the  experimental  effort. 

1.  Design  of  the  test  rig  has  been  completed,  and  various  dimensions 
and  materials  for  it  have  been  selected,  guided  by  analysis  of 
various  geometric  inputs  on  gas  dynamic  operating  condition 

in  the  test  rig. 

2.  Fabrication  of  major  components  of  the  test  rig  has  been  completed. 

3.  Various  instruments  and  other  accessories  have  been  selected  and 
purchased. 

4.  Propellant  grain  for  the  driving  motor,  test-propellant  strips, 
and  pyrotechnic  igniter  have  been  acquired. 

5.  A four  channel  electronic  control  unit  for  the  stepper  motors 
has  been  designed  and  fabricated. 

6.  Some  tests  of  the  laser-photodiode  servomechanism  were  conducted. 
Feeding  mechanism  and  electronic  integrading  units  operated  properly, 
in  order  to  record  feeding  rate  of  burning  propellant  strips. 

These  same  preliminary  tests  illustrated  the  feasibility  of 
filtering  light  emission  from  burning  propellant  via  narrow  band 
pass  filters.  Laser  beams  can  therefore  be  isolated  and  detected 
by  photodiodes. 

IV . C . Flexibility  of  the  Experimental  Apparatus 

The  experimental  apparatus  for  erosive  burning  studies  has  the 
following  flexibilities: 

1.  Free  stream  velocity  and  mass  flux  through  the  test  chamber 
can  be  varied  by  changing  the  exit  nozzle  throat  area  and  the 
height  of  the  top  plate  in  the  test  chamber. 

2.  Pressure  in  the  test  chamber  can  be  controlled  by  using  different 
exit  nozzles. 
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3.  Pressure  gradients  in  the  streamwise  direction  can  be 


controlled  by  using  a tapered  interchangeable  top  plate. 
4.  Propellant  strips  of  different  compositions  can  be 
tested. 


42 


V.  WORK  TO  BE  PERFORMED 

Theoretical  Work 

The  following  work  is  to  be  done  in  this  portion  of  the  current 
research  program: 

1.  Turbulence  Modeling:  Modeling  of  the  turbulence  correlations  listed 

under  section  III.C.  is  in  progress.  These  correlations  will 

be  related  in  terms  of  known  variables.  Turbulence  transport 
phenomena  and  the  gas-phase  combustion  will  be  linked  through  modeling 
of  the  correlations.  The  boundary  conditions  required  for  the 
turbulence  related  terms  ( k and  UV'  ) are  also  in  progress. 
Introduction  of  particle  size  effect  on  erosive  burning  through 
the  boundary  conditions  on  U.V'  is  under  consideration. 

2.  Numerical  Solution  and  Parametric  Studies:  As  soon  as  the 

theoretical  formulation  has  been  completed,  three  major  tasks 
in  the  numerical  solution  will  begin.  The  first  task  is 
selection  of  a stable  and  accurate  numerical  scheme.  Literature 
survey3 0-3 3 is  already  in  progress  to  accomplish  this  (Crank- 
Nicholson  implicit  scheme30  is  under  consideration).  The  second 
task  is  computer  program  implementation  and  debugging.  To 
facilitate  debugging,  coding  will  emphasize  logical  simplicity 
rather  than  computational  efficiency.  Computational  efficiency 
can  be  improved  readily  once  program  validity  is  established.  The 
third  task  is  the  convergence  tests  and  the  error  estimations. 
Convergence  tests  will  be  carried  out  by  varying  the  mesh  size. 

The  solutions  must  change  very  slightly  when  the  mesh  sizes  are 
reduced.  Since  convergence  of  numerical  solutions  depends  upon 
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how  well  the  problem  is  posed,  various  tests  can  be  performed 
on  solution  sensitivity  to  changes  in  boundary  conditions,  and 
to  changes  in  other  important  parameters  such  as  mass  flow 
velocity,  chamber  pressure,  free  stream  temperature,  initial 
temperature  of  the  propellant,  oxidizer  particle  size,  oxidizer 
to  fuel-binder  ratios,  and  strand  burning  rate  parameters  of  the 
propellant.  In  all  these  tests,  if  the  solutions  are  bounded  and 
change  only  slightly  for  small  perturbations  on  the  parameters, 
the  existence  of  the  neighboring  solutions  is  proved,  and 
convergence  of  numerical  solutions  is  assured.  Furthermore, 
convergence  and  accuracy  of  solutions  can  be  verified  by  carrying 
out  residual  iterations,  and  the  magnitude  of  the  truncation  errors 
can  be  determined. 

Once  the  numerical  evaluation  technique  is  finalized,  a 
series  of  "computer  experiments"  will  be  performed  to  determine  the 
influence  of  some  particular  parameters  of  interest  such  as 
combustion  gas  velocity,  flame  temperature,  pressure  level, 
weight  fraction  of  fuel  binder,  oxidizer  particle  size,  etc. 

The  results  will  be  used  in  generating  an  expression  for  burning 
rate . 

Experimental  Work 

As  soon  as  the  fabrication  of  the  test  rig  has  been  completed,  a 
series  of  experiments  will  be  conducted.  The  experimental  measurements 
will  include: 

1.  burning  rates  of  propellant  strips, 

2.  mass  flow  rate  of  the  combustion  gases  in  the  test  chamber. 
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3-  pressure  and  pressure  gradient  in  test  chamber  and, 

4.  temperature  of  the  combustion  gases. 

A high  speed  movie  camera  will  be  used  to  observe  the  burning  process 
during  an  experimental  firing. 

Development  of  an  Erosive  Burning  Rate  Equation 

Controlling  dimensionless  variables  will  be  identified  from  the 
theoretical  model  and  dimensional  analysis.  The  variables  will  be  measured 
directly  or  deduced  indirectly  from  the  experiments  and  effect  of  each  on 
erosive  burning  will  be  studied  experimentally.  This  information  will  be 
used  to  determine  the  functional  form  relating  all  the  dimensionless 
parameters.  Numerically,  parametric  study  will  be  carried  out  to  extend 
the  range  of  variation  of  these  parameters. 

Measured  and  numerically  calculated  values  of  these  parameters  will 
be  fed  into  a multiple-regression  analysis  program  to  determine  the  final 
correlation  for  erosive  burning  rate  in  terms  of  the  non-dimensional 
parameters.  The  multiple  correlation  coefficient  and  standard  deviations 
of  the  exponents  and  constants  will  be  evaluated, in  order  to  indicate  the 
degree  of  correlation  with  experimental  data.  The  erosive  burning 
expression  will  be  useful  for  propellant  grain  designers,  who  wish  to 
achieve  better  control  of  erosive  burning  behavior  in  high  performance 
rocket  motors. 
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APPENDIX  A 


ORDER  OF  MAGNITUDE  ANALYSIS 


For  the  order  of  magnitude  analysis  of  Eqs.  (8)  - (12),  the  dependable 
variables  will  be  referred  to  their  value  at  some  reference  point,  \ , 
outside  the  boundary  layer,  and  coordinates  x and  ^ will  be  referred, 
respectively,  by  a representative  dimension.  Thus,  we  write 

U - 4 0; 

1 X 

■ax  L bf  S 


Order  of  magnitude  of  p'u.'  and  f7?  : The  density  fluctuations,  y'  , 

result  from  several  causes26’27  listed  below. 

(a)  Mixing  of  fluid  elements  by  turbulent  bulk  transport  from 
regions  of  different  density,  velocity  and  enthalpy. 

(b)  Turbulent  pressure  fluctuations — compressibility  effect. 

(c)  Molecular  viscous  dissipation  caused  by  fluctuating  velocity 
and  molecular  heat  transfer  driven  by  the  fluctuating  enthalpy. 

According  to  references  23,24,26,27,  (a)  is  considered  to  be  the  dominant 

cause  for  production  of  f'  . We,  therefore,  assume  f7  — l,  , and 

0 

similarly  where  l, , (t  are  turbulent  length  scales.  Assuming 

l,  and  £z  of  the  same  order,  we  get 

ps.f7 . £>(  SV-§-) 

n 


= 0(C7 


aF 

■Jr) 

6 
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W * f'v'  = 0 


(A.  1) 


From  the  x-momentum  equation  (Eq.  9),  If  turbulent  shear  stress  is 
important , we  must  have 

UW  - o(fE|s) 

But  (fvVU'  ^ 0 ( F ^ ) 

7,0^)]-  o(?ttl&) 

U^jOt1  + °(?  lO 

In  the  present  problem,  temperature  ratio  across  the  boundary  layer  is  of 
the  order  of  2 to  3,  therefore  _ 

4^  oo) 

X 

(a. 2) 

Substituting  Eq.  (A. 2)  into  Eq.  (A.l),  we  get 

fuT~  pV-0  (A- 3) 

Order  of  magnitude  of  f'h  : Kistler2 1 has  found  experimentally  that  the 

distribution  of  static  temperature  fluctuations  is  very  nearly  proportional 
to  the  static  temperature  difference  across  the  boundary  layer.  Therefore, 
we  can  assume  that  the  enthalpy  fluctuations,  ^ , result  mostly  from  the 
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turbulent  bulk  transport  of  fluid  elements;  we  can  then  write 


Assuming  the  length  scales  l, , tt  and  li  are  of  the  same  order  of  magnitude, 
we  get 


/ 


iL  dL-  x 


, . d at  V 

s1 


(A. 4) 


(A.  5) 
(A.  6) 

In  the  following,  the  steady  state,  two-dimensional  conservation  equations 
are  written  and  the  relative  orders  of  magnitude  of  the  terms  are  given 
below  each  equation. 

Mass  Conservation  Equation: 


Similarly,  we  can  show  that 

^ o (i  a*17*) 
0(£ 


o * 

- k<p7>  + vr,) 

f = o 

(A. 7) 

C‘*f> 

(ff)  0-f) 

cf) 

x-Momentum 

Conservation  Equation: 

f>a$  - 

• - 

r?lf 

0) 

(*£  f)  0) 

ff) 

ii 

i 

+ * [J-  % * i is, 

+.  3j? 1 — ?_Cfu)'u,' 

T>* 

- |j(WV  (*.8) 

M?  »*' 

O («♦*>* 
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I 

I 


y-Momentum  Conservation  Equation: 


faff  + 

f’S'g  * - 

7752 

(4) 

(i) 

(f  tt 

W 

i 

ii 

+ /*■  [ IS  "•  3 tftj  + 

i i2]  . 

_ 2_  (fu)V 

_ 5-  cfy/v' 

L 

(’•  pi’ 

C*f) 

Energy  Conservation  Equation: 

After  replacing  the  viscous  dissipation 
dominant  part,  />-  (5^)  » E9*  (11)  becomes 

term, 

, by  its 

fag 

+ pn?||  * 

+ 

77  aK 

a) 

(f{)  (0 

<f> 

_U  *§ 
ax 

rt  al»  a>  ( > &K 

" v ay  c»>  a* 

) + 

5.  ( 1 ?L-\ 

- 

/ L <>*n 
D*  / 

{L  L ( 1 >1 

L j>  a|j  1 1 R.P^ 

RttW 

r«  + 

L T 

- lu(p/)^ 

O.f) 

,Lf  1 B? 

's'  5 if 

I 

I 

I 


(A. 9) 


(A. 10) 


Reynold's  Stress  Equation: 

The  terms  u.'} -*■  of  Eq.  (12)  can  be  simplified  as 

follows 

Let  [I]  = Ujrlj.j  +uT^j.( 

using  T'j  = M ( -t-  Uj,i  ^ U'*r,-m') 

we  find 

[I]  = Uj  CMUlj-i.j  + UiC/UU',,p^  + 

— LL'iCM  LLim, ■***),  I 

Let  us  assume  (for  finding  relative  importance  of  various  terms  in  I ) 
that  /Um  constant  in  the  last  four  terms.  Therefore,  we  can  write 

[I]  * ”*/*  UW“'7  +J  + J (A- 12) 

Next,  we  will  show  that  the  last  two  terms  of  Eq.  (A. 12)  are  small 
and  can  be  neglected.  To  find  the  order  of  magnitude  of  these  two  terms, 
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differentiate  the  mass  conservation  equation  of  fluctuating  part  with 
respect  to  xi  , multiply  the  result  with  u£  , and  finally  average  the  whole 
equation.  The  result  is  given  by  the  following  equation: 


Ui  I-.,# 


+ f U; 


+ + u-i  c f — o 

At  the  most  p uTUj.ji  ^ 0 ( ) 

-°(fP 

Here  «-Tis  a characteristic  turbulence  velocity. 
Now  in  Eq.  (A.  12)  Ui.j  u'f.j  ^ Q ( &) 

where  ^Tis  the  Taylor  microscale 


/2rf  el 


« 


because  ^t<<8 


we,  therefore. 


can  write 


^ SUf 


(A. 13) 


With  this  simplification  and  the  assumption  of  steady  state,  we  can  write 


Eq.  (12)  into  the  following  two  equations: 

Turbulent  Kinetic  Energy  Equation:  0=0 


(fu = 


- (fup'ui  Ui,j  + 


- <7^-  - O' 


(A. 14) 
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Reynold's  Stress  Equation:  (x=i,  f=i) 


ffOj+  f! j)  ^ Wj  - «’“PV 


* + <•'<?$'  - £)  - 


+ 


- 2M 


au'  W 

*V'i 


(A. 15) 


Procedure  for  the  order  of  magnitude  analysis  of  Eqs.  (A. 14)  and 
(A. 15)  is  similar  to  that  discussed  earlier.  The  result  is  that  the 
terms  with  mean  gradient  quantities  in  x-direction  are  negligible. 

After  a comparison  is  made  between  the  various  orders  of  magnitude 
in  each  equation,  Eq.  (A. 7)  - (A. 11),  Eq.  (A. 14)  - (A. 15),  only  the  dominant 

terms  are  retained.  It  is  noted  that 

<:<  1 > Re  > 7 i 

(—■f  $ i » ft*  1 f Se  ~ 1 

Finally,  we  obtain  the  following  boundary  layer  equations  for  a 

chemically  reactive,  compressible,  turbulent  boundary  layer  developed 

over  a flat  plate: 

Continuity  Equation: 

fjC?a)  + V”)  = o 

where  s fV  (A. 16) 

Momentum  Equation: 
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Species  Equation: 
Energy  Equation: 


pall?  - ?''r»  = * a« 


+ |j[-  (WF  +W]  - + T^U \ti 

Turbulent  Kinetic  Energy  Equation: 

fag  f5<|k  = + P'Z.i 

-|»Cf<k'*f5^  - /“  SSJ TZ,i 

Reynolds'  Shear  Stress  Equation: 


f U + f\P|-(W)  = 


- tjcfcA  f'KJ  * 

-ZMWW 


<)xj 


(A. 18) 


(A. 19) 


A. 20) 


(A. 21) 
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APPENDIX  B 


BASIC  EQUATIONS  AND  PARAMETRIC  STUDIES  IN 
THE  TEST  RIG  DESIGN 


The  test  rig  consists  of  four  major  components:  the  driving  motor, 

the  test  chamber,  the  feed  mechanism  chamber  and  the  pressure  equalizer 
chamber  (see  Fig.  B.l).  The  design  of  these  components  of  the  test  rig 
depends  on  the  gas  dynamic  operating  conditions  for  pressure,  temperature 
and  velocity  required  in  the  test  chamber.  Therefore,  the  equations 
of  mass  and  energy  balance  in  the  driving  motor  and  the  test  chamber 
are  considered.  Through  these  equations,  the  effect  of  various  geometric 
inputs  of  the  test  rig  on  the  gas  dynamic  operating  condition  can  be 
studied. 

The  equations  for  mass  and  energy  balance  in  the  driving  motor  and 
the  test  chamber  can  be  derived  by  considering  a control  volume  in  each 
of  these  components,  and  introducing  the  following  assumptions: 

1.  All  the  properties  are  functions  of  time  only,  and 

2.  Perfect  gas  law  is  valid. 

Mass  balance  in  the  driving  motor: 

Energy  balance  in  the  driving  motor: 

I +"  ('VT,'»c4>Tl)d.  “ -3)4  — 


(B.l) 


(B.  2) 


where  the  subscript  d refers  to  the  driving  motor, 
using  M = P V 

and  perfect  gas  law  b = -^TT‘  ’ e = ‘ 

w 

we  find, 

dM  n dV  . H ..  ±1  4Z 

it”  ' it  RT  dt  T dt 


(B.3) 


(B.4) 
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From  Eqs.  (B.3)  and 
d(Mf)  _ r PT 


Substitute  Eqs.  (B.4)  and 
have 


(B.4),  we  have 

dV  CvVd+> 

It  + R dt 

(B.5)  into  Eqs 


(B.l)  and  (B.2) 


respectively,  we 


( ? 51  + & ft  ~ ^ + <*<»>*  -(■*«#>* 

(CvfTg  - ^ $ ld  = 


In  the  driving  motor,  the  rate  of  change  of  volume  is  given  by  the 
following  equation: 

(£)„= 

Multiply  Eq.  (B.6)  by  (CyT)^  and  subtract  from  Eq.  (B.7),  and  use 
Eq.  (B.8),  we  finally  get 

/dl\  = 

-Wi  - (S-+ 


Putjr)  -v  C^b  — C^exit 


C™b^Tf^-(n"W  ^TfxJdi  — 


(B.5) 

(B.6) 

(B.7) 

(B.8) 

(B.  9) 
(B. 10) 
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Mass  balance  in  the  test  chamber: 


(^r)c=  + C™b\  ~ 0*st,i«Oe 

Energy  balance  In  the  test  chamber: 

(ir^lc=  -K^b CV  Tf  )t  - (tvw  C™s**f 


The  subscript  c refers  to  the  test  chamber.  A simplification  of  Eqs.  (B.ll) 
and  (B.12)  similar  to  that  followed  earlier  for  the  driving  motor  results 
in  the  following  equations: 


(dfX-  (“ipRV")c  [ +C'roh  ^ ^ 


-CW>C-  | (tt)c^"c 


(gt)c=  |c  Verify  tn-s,,,Ot 

-cU  - (gr)J%\T*  ] 


In  Eqs.  (B . 9)  - (B.10)  and  Eqs,  (B.13)  - (B.14),  <vn(!xif  . and  4vyexj- 
are  calculated  by  the  following  equations: 


(B.ll) 

(B.12) 


(B.13) 

(B.14) 


(B.15) 
(B. 16) 
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Eq.  (B.15)  is  based  on  Che  assumption  of  isentropic  flow  in  a 
nozzle . 

We  can  also  write  for  the  volume  in  the  test  chamber, 

Vc  = A* (B.  17) 

& = (B.  18) 

where  (Vo)c  is  the  volume  of  the  test  chamber  when  x^,  =.  o (see  Fig.  B.l). 

and  depend  on  the  pressure  difference  across  the  piston.  The  pressure 
equalizer  chamber  consists  of  2 cylinders.  Each  cylinder  contains  a 
piston,  a spring  and  a damper.  A force  balance  for  the  piston,  spring 
and  damper  system  within  a single  cylinder  gives  the  following  equations: 


— {”(£  A*  - dlcU*  ~ (B.  19) 

Oil  Vr\p 

= U-4>  (B.  20) 


where  is  the  piston  mass,,  c*c  is  the  damping  factor,  dv.  is  the  spring 
constant.  To  evaluate  , we  assume  a poly tropic  compression  of  the 
gas  in  the  feed  mechanism  chamber.  Therefore,  we  have 

<R»o.  T (8.2i) 

where  (Vo)^  is  the  volume  in  the  feed  mechanism  chamber  excluding  the 
free  volume  in  cylinders,  is  the  initial  pressure,  -r\  is  the 

exponent  in  the  polytropic  compression, and  r»c  is  the  number  of  cylinders. 

Surface  temperature  of  propellant  strips  in  the  test  chamber  is  given 
by  the  following  equation.* 


cdt  '*■  3^(WUii-vV) 


*This  equation  is  derived  from  the  transient,  one-dimensional  heat  conduction 
equation  in  solid  propellant  (see  reference  37). 
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where  ^v0»>v  is  given  by  the  Dittus  Boelter  relation: 

*Wv=  Pv5 


"ffj  is  assumed  constant  after  it  has  equalled  the  ignition  temperature. 

Eqs . (B. 8)  - (B.10),  (B.13)  - (B.14),  and  (B.18)  - (B.22)  have  been 
solved  on  the  computer  using  fourth  order  Hamming's  Modified  Predictor- 
Corrector  Scheme.  The  properties  used  in  the  calculations  are  listed  in 
Table  B.I.  A parametric  study  has  been  performed  to  determine  the  effect 
of  various  geometric  inputs  on  the  gas  dynamic  properties:  pressure, 

temperature,  velocity,  etc.,  in  the  test  rig.  The  results  of  these 
calculations  have  been  plotted  and  shown  in  Figs.  B.2  - B.7. 

Conclusions  drawn  from  the  parametric  study  are: 

1.  Various  gas  velocities  (0-250  m/sec)  and  various  pressures  (up 
to  300  atm)  can  be  achieved  in  the  test  chamber  by  controlling 
exit  nozzle  area  (see  Figs.  B.4  and  B.5). 

2.  Velocity  can  also  be  controlled  by  the  variation  of  height  in 
test  chamber  (see  Fig.  B.6). 

3.  Mass  flux  and  velocity  in  the  test  chamber  can  be  changed  by 
using  different  burning  surface  area  of  propellant  in  the  driving 
motor  (see  Fig.  B.5  and  B.7). 

4.  Steady  state  conditions  for  pressure,  temperature,  and  velocity 
can  be  achieved  in  the  test  chamber  after  a very  short  transient 
interval  (see  Figs.  B.2  and  B.3). 
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TABLE  B.I . 


Properties  Used  In  Parametric  Studies 


(W)d  - 22.01  Kg/Kg-mole 


(W)  - 22.01  Kg/Kg-mole 

c 


<P.>d 


1.6  x 10~3  Kg/cm3 


(p  ) - 1.6  x 10"3  Kg/cm3 

s c 

(as  ) - 1.875  x 10  3 cm^/sec 

C 

(As  )e  ■ 9.0  x 10  7 Kcal/cm-0K-sec 
-ft  —TO 

U - (2.5  x 10  + 4.333  x 10  T ) Kg/cm-sec 

c c 

(T . ) - 700°K 

ign  c 
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Fig.  B.3  Calculated  Temperature-Time  Traces  in  the  Test  Rig 


Pressure  In  Test  Chamber 


Average  Velocity  in  Test  Chamber 


n ‘aaqnreqo  isax  uf  ^3700x3^  33bj3av 


Fig.  B.5  Effect  of  Burning  Surface  Area  in  Driving  Motor  and 
Exit  Area  on  the  Velocity  and  Pressure  in  the  Test 
Chamber 
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Fig,  B.7  Effect  of  Burning  Area  in  Driving  Motor  and 
Exit  Area  on  the  Mass  Flux  and  Velocity 
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Velocity  in  Test  Chamber,  U , m/sec 
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